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A)  OBJECTIVE. 

Determination  of  the  location  on  the  earth’s  surface  of  an  electromagnetic  emitter  by  a 
single  satellite  sensor  or  a  network  of  satellite  sensors,  is  possible  if  the  variables  which  define  the 
location(s)  and  orientation(s)  of  the  sensor(s)  is  known  exactly.  However  if  an  emitter  location 
is  determined  according  to  inaccurate  information  of  sensor(s)  location  and  orientation,  then  the 
determination  will  also  be  in  error. 

In  this  study,  we  will  first  handle  the  case  of  a  single  sensor.  Networks  of  two  or  more 
sensors  will  be  handled  later.  The  error  in  the  location  determination  is  defined  as  follows:  In 
figure  1,  let  5,  represent  the  location  and  orientation  of  the  sensor  according  to  inaccurate 
information.  Also  let  Ix  represent  the  emitter  location  due  to  that  data.  Let  S2,and  I2  be 
identified  as  were  Sv  /,  but  according  to  exact  data.  The  error  is  then  the  vector  difference 
I2-Ir  In  this  study,  I2  -Ix  will  be  approximated  by  its  total  differential  dl .  The  formula  relating 
the  total  differential  to  the  sensor  variables  will  be  used  to  find  the  absolute  maximum(s)  of 
\dl\  with  respect  to  changes  in  the  sensor  variables. 
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Figure  1.  Geometry  of  the  relation  between  the  error  I2  -Ix  and  the  assumed  and  true  satellite 
positions  5,  and  S2  respectively. 
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Definition  of  the  sensor  location  and  orientation  variables 


Referring  to  figure  2,  consider  the  following  situation.  A.satellite  sensor  system  is  located  at 
point  S  at  a  certain  instant  of  time  tQ .  The  sensor  is  orientated  so  that  at  tQ ,  it’s  “sight  axis” 
is  along  SI  which  meets  the  earth’s  surface  at  I.  At  that  time  the  sensor  detects  a  signal  from  I. 


Figure  2.  Geometric  relationship  between  the  sensor  at  S  and  an  emitter  at  I  on  the  earth’s 
surface. 
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Introduce  the  inertial  geocentric  right  hand  coordinate  system  [Gj  ]  with  plane  i(G  ^  -  y[G^ 
coincident  with  the  equatorial  plane."  It  will  be  assumed- that  system  [Gj]  is  the  known  basic 
coordinate  system  into  which  all  vectors  will  be  resolved  in  order  to  compare  to  other  vectors. 

The  vector  GS  from  geocenter  to  satellite  is  described  in  terms  of  the  three  quantities, 

A ,  P  and  r  -  re+h  which  are  measured  as  indicated  in  figure  2.  The  angles  ap  and  Pp  then 
complete  the  description  of  the  sensor’s  attitude.  These  last  two  angles  define  the  orientation 
with  respect  to  system  [Sj] ,  which  is  characterized  by  having  origin  at  S,  axis  ks^  colinear  with 
GS  and  axis  pointing  North  (see  figure  2).  Thus  the  five  variables  r,  A,  P,  ap,  and  Pp  are 
sufficient  to  describe  the  location  and  orientation  of  the  sensor. 


■a  System  [G,]  has  coordinate  axes  i[C^,  j[G^ ,  k[G  ].  This  convention  is 
used  throughout  this  paper  so  that,  for  example  coordinate  system  [LJ  has  coordinate  axes 
hLt]  »  j[Lkr  ^  coordinate  systems  will  be  right  hand  cartesian.  Furthermore  the  title  of  any 
particular  coordinate  system  will  indicate  its  origin  of  coordinates  so  that  for  example,  all 
[Gj\  systems  will  be  geocentric,  while  alltSJ  systems  have  origin  at  S.  Also,  letters  in 
brackets  will  indicate  coordinate  systems.  Finally,  a  quantity  that  has  a  subscript  consisting  of  a 
bracketed  symbol,  will  indicate  the  representation  of  that  quantity  in  that  coordinate  system.  Thus 
e.g.  I[Si]  is  the  representation  of  I  in  [Sj]  coordinates. 


4 


ANALYSIS. 


The  detailed  geometric  action  of  the  sensor  can  be  thought  of  as  a  mapping  of  a  point  I  on 
the  earth’s  surface  onto  an  image  point  on  the  z-axis  of  a  coordinate  system  with  origin  at  S  as 
follows:  Refer  to  figure  3  and  introduce  coordinate  system  [S2]  which  is  characterized  by  having 
axis  opposite  to  SI  and  axis  j[S^  parallel  to  SPxSI 


Figure  3.  Geometry  of  Coordinate  System  [S2] 
.a  SP  is  the  projection  of  SI  on  the  plane  i^,  j[S^ 
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Now,  knowing  the  quantities  r,  X,  p,  a p,  Pp,  we  can  determine  the  location  of  S  and  the 

#  1 

representation  of  vectors  SI  and  SP  in  [GJ  (See  figure  2).  Then  as  these  last  two  vectors 
determine  two  axes  of  coordinate  system  [S2] ,  it  is  possible  to  obtain  completely  the  rotation 
matrix  (M)  from  [S2]  to  [SG]a.  The  complete  transformation  of  point  “I”is  then  the  result  of  a 
translation  along  IS,  a  rotation  through  the  matrix  (M)=(m~)  where  m- -f^iX,  p,  dp,  P^)  anda 
translation  along  SG.  If  now  the  quantities  r,  X,  p,  pp  are  altered  by  their  respective 
increments  to  values  r',  X',  p7,  dp,  p^,  (where  for  example,  r'=r +Ar),  S  shifts  to  S'  ,  and  the 
vector  SI  becomes  S  'l ’  (see  figure  4).  The  error  l'-I  will  be  approximated  by  the  total 

differential  dl  represented  in  coordinate  system  [G,] . 


■a  The  system  [SG]  is  parallel  to  system  [GJ  ,but  with  origin  at  S. 
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The  process  of  determining  the  point  I  on  the  earth  and  also  its  total  differential  dl  will  now  be 
described  in  detail.  With  S  as  the  origin  of  the  [S2]  coordinates,  the  [S2]  coordinates  of  S  are 
of  course  (0,0,0 ).  We  first  determine  the  point  I,  which  is  the  intersection  of  SI  with  the  earth  s 
surface.  In  [S2]  coordinates,  the  coordinates  of  I  are  (0,0,1  Ik^  is  determined  in  Appendix 
II.  Let  (M)  be  the  rotation  matrix  that  maps  coordinates  in  [S2]  to  those  in[SG]  (introduced 
previously).  (M)  is  derived  in  appendix  1.  Thus  the  [SG]  coordinates  of  I  are 
/  =  (M) (0,0,1.  )T.  Finally,  a  translation  by  re  +  h  along  the  negative  of  yields  I  in 

[G,]  coordinates.  Finally,  we  compute  dl,  the  total  differential  of  I  with  respect  to  the  variables 
that  define  sensor  location  and  orientation,  dl  is  the  representation  of  the  error  that  we  seek. 
Thus,  given  a  set  of  values  of 

r,  X,  p,  ap,  pp,  then  the  location  of  the  emitter,  I  is  =  Oi^  +  0 ;[5J  - 
where 

\SI\  =  rsinPp-^/-e2-(rcospp2 

as  shown  in  Appendix  H.  Then  pre-multiplying  Z[5J  by  the  rotation  matrix  (M)  (derived  in 
Appendix  I)  yields  I[Sc] 

bd  =  W  hs2] 

as 

I[SC]  =  I  SI\(mni[Sci  +  m23 j[Sc]  + 

Now  doing  a  translation  along  the  vector  GS  (which  is  along  the  fc[S;]  axis)  and  using  the  first 
equation  in  Appendix  I,  we  obtain 


7 


ho,]  =  hsci  - 1 SG I  (i[5c]cosPcosA  V^cos^smX  +^cJsinp) 
which  upon  substitution  from  the  above  can  also  be  written 

7[Cj]  =  (m13|S/|-|SG|cospcosX)i[0)  ♦  (m23  |S/|  -  |SG|  cosp  sinA)j(G1  , 

(^IS/HSGIsinP)*^ 

where  in  the  immediately  above  relation  we  have  recognized  that  the  coordinates  systems 
[Gj]  and  [SJ  have  identical  coordinate  axes.  This  constitutes  the  desired  quantity  and  coordinate 
system  in  which  to  represent  it.  Next,  the  total  differential  dl  with  respect  to  the  variables 
r,  X,  P,  a p,  pp  and  the  function  |  dl  \ 2  will  be  formed  The  first  partial  derivatives  of  \dl\2 
will  be  formed  with  respect  to  the  variables  r,  X,  p,  a^,  Pp  and  critical  points  identified.  Domain 
end-points  will  also  be  inspected.  From  these  points,  absolute  maximum(s)  of  |  dl\2  will  be 
determined. 
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APPENDIX  1 


DERIVATION  OF  THE  (M)  MATRIX 


In  order  to  derive  the  (M)  matrix  it  is  necessary  in  what  follows  to  work  with  coordinate 
systems  having  a  common  origin.  Accordingly  we  have  previously  introduced  the  system  [SG]  as 
being  parallel  to  [G,  ]  but  with  origin  at  point  S. 

By  recalling  Figure  2  and  our  definition  of  system  [Sj  ] ,  it  is  clear  that  the  representation  in 
coordinate  system  [SG]  of  axis  is 

\s,)  =  hs0]c ospcosA  +  y[5c]cosPsinA  +  A^sinp. 

Furthermore,  axis  i[s  ]  points  North  from  point  S  and  therefore  lies  in  the  plane  containing 
axes  k[5  ]  and  ]c[G ^ .  Then  as  axis  j[S  ]  is  perpendicular  to/[S  ]  and  also  system  [5,]  is  right-hand, 


we  have  for  the  situation  of  Figure  2 

fc[S,]x  km 

[5/]X  klG,l 


=  r[5c]sinX  +j[S/~  cosX)  +  ^0 


Finally,  in  completing  the  right-hand  system  [SJ  ,we  construct 


hs,)  =  =  ilSGi(-cos\sm  P)  +  ./^(sinAsinP)  +  ^cosp. 

Then  the  matrix  that  yields  the  components  of  a  vector  represented  in  system  [SG]  in  terms  of  its 
components  in  system  [5J  is 

-cosAsinP  sinA  cosPcosA^ 

(o,)=  -sinAsinP  -cosA  cosPsinA 

>  cosp  0  sinp 
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The  remainder  of  the  derivation  proceeds  in  an  entirely  analogous  way  from  the  definition  of 
coordinate  system  [S2]  described  in  the  opening  paragraph  of  the  section  labeled  “analysis”. 

Let  lSI  be  a  unit  vector  in  the  direction  of  SI.  Then  from  Figure  2,  we  have  in  [SJ  coordinates 

*$/  =  i[s,)cos  P/>cosap  +  j[Sj]cos$psmap  +  £[5J(-sinpp). 

By  definition,  the  representation  of  axis  k[S^  in  system  [SJ  is 

kis2]  =  ~hi  =  ^(-cosp^cosa^)  +7[5](-cosppsina/))  +  £[5/]sinpp 

Furthermore,  is  parallel  to  SPxSI  and,  with  l5p  represented  in  system  [5,]  (See  Figure  2) 
as 

=  l[S7]COS  aP  +  j[S1]S^aaP  +  ^ISj]^ 

we  obtain 

•  I5/ 

hs2]  =  J^x  t ^  =  hs/~sinap)  =%( cosap  +  k[S]] 0 
Finally,  since  [S2]  is  to  be  right-hand, 

hs2]  =  j[s2]x  kis2]  =  Ic5;]cosapsinPp  +j\sj]  sina^sinP^  +fc[5;]cospp 
Thus  the  matrix  (a2)  from  system  [S2]  to  system  [5.]  is 

I  1 

/CosapsinPp  -sinap  -cosapcospA 

(o2)  =1  sinapsinpp  cosap  -sinapcospp  I 


so  that  the  matrix  (M)  from  [S2]  to  [SG]  has  the  form  (M)  =  (a,)(o2)  With  the  following 
elements: 

mn  =  -  cos  A  sin  P  cos  o^sin  Pp + sin  A  sin  a^sin  Pp + cos  p  cosA  cos  P p 
mn  =  cosXsinpsinap+sinXcosap 

mn  =  cosAsinPcosa/>cosP/>-sinAsina/,cosPf>+cosPcosAsinP/> 
m2l  =  -  sin  A  sin  p  cos  apsin  Pp  -  cos  A  sin  a^sin  Pp + cos  P  sinA  cos  pp 
m22  =  sinXsinPsina^-cosXcosap 

m23  =  sin  X  sin  P  cos  apcos  $p + cos  X  sin  a^cos  P^ + cos  P  sinX  sin  pf 

m31  =  cos  P  sin  p^cos  ap  +  sin  p  cos  Pp 

m32  =  -cos  P  sin  ap 

m33  =  -  cos  P  cos  PpCos  ap + sin  P  sin  Pp 
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In  order  to  do  this,  we  first  represent  the  vector  GI  (from  geocenter  to  point  I)  in  [S2] 
coordinates  as  G/^  ]  =  GS^  +  Sfs2]  In  system  [S2] ,  GS  has  the  representation 

GS[s2]  =  I[sjrcos  +  j[s2l®  +  hsf^p 

also  we  have 

G/[SJ  =  ift/00*  +  As,)0  *  V“  lS/l  +rsinP/'> 

(With  the  minus  sign  in  front  of  |  SI  |  because  SI  is  along  the  negative  of  the  axis)  with  the 

value  of  |S/|  to  be  determined. 

As  I  is  on  the  earth’s  surface,  \GI\  =  the  radius  of  the  earth  re.  Inserting  this  condition  and 
forming  |  GI\ 2  we  obtain  the  quadratic 

\SI\2  -  2|S/|rsinPp  +  r2  -  r2  =  0 

with  the  solution 

\SI\  =  rsinpp  -  Jr2  -  (rcospp)2 

where  the  negative  square  root  has  been  chosen  because  it  corresponds  to  the  smaller  value  of 
\SI\  ( since  0<  Pp<  7t/2).  Figure  5 reveals  that  the  positive  square  root  corresponds  to 
extending  line  SI  through  the  earth  and  intersecting  it  on  the  other  side. 
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SUMMARY  OF  THIS  REPORT  AND  OUTLINE  FOR  WORK  STILL  TO  BE  DONE 


The  development  of  an  equation  relating  the  location ‘T’  of  an  electromagnetic  emitter  to  a  set 
of  variables  defining  location  and  orientation  of  a  single  satellite  sensor  system,  has  been 
accomplished.  Next,  the  total  differential,  dl  with  respect  to  the  variables  r,  X,  p,  P p,  ap  will  be 
formed  as  will  the  function  \dl\2.  For  this  last  function,  the  absolute  maximums  will  be 
determined  by  examining  critical  points  and  domain  end-points.  This  constitutes  the  work  for  the 
single  sensor  system. 

With  the  pattern  of  analysis  defined  by  the  single-sensor  system  above,  the  case  of  a  two-sensor 
.  system  will  add  five  more  variables  rv  X2,  P2,a2^,  P2^and  a  filtering  scheme  with  which  to  combine 
the  output  of  the  two  sensors.  Three  or  more  sensors  are  handled  in  an  exactly  analogous  manner 
as  was  the  two  sensor  system  described  above. 
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